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ABSTRACT 



This study presents the derivation of a diffuser design parameter, 
which takes into account losses due to the geometry of the diffuser. 

A criterion for design of annular diffusers, using the design 
parameter to predict losses, is developed from reference data and 
model tests . 

A curved wall annular diffuser was designed and tested to verify 
the predictions of the design criterion for an arbitrary diffuser. 



2 



TABLE OF CONTENTS 



I. INTRODUCTION 13 

II. DEVELOPMENT OF THE SHAPE PARAMETER OMEGA (fi) 16 

A. ENTROPY CHANGE OF A ONE -DIMENSIONAL, FRICTIONAL FLOW . . 16 

B. ENTROPY CHANGE FROM THERMODYNAMIC CONSIDERATIONS .... 19 

C. DEFINITION OF OMEGA (Q) 21 

D. REIATIONS BETWEEN Q AND DIFFUSER PERFORMANCE 22 

III. TEST INSTALIATION AND INSTRUMENTATION 26 

A. AIR DELIVERY SYSTEM 26 

B. DESCRIPTION OF DIFFUSER MODELS AND ASSEMBLY 28 

C. INSTRUMENTATION 33 

IV. COMPUTATIONS AND EXPERIMENTAL PROCEDURES 38 

A. PRESSURE 38 

B. TEMPERATURE 38 

C. FLOW RATE 39 

D. CHECK OF NOZZLE FLOW RATE WITH SURVEY DATA 40 

E. REYNOLDS NUMBER 44 

V. ANALYSIS OF TEST DATA FROM REFERENCE 2 45 

VI. DISCUSSION OF TEST DATA 48 

VII. CONCLUSIONS AND RECOMMENDATIONS 60 

APPENDIX A - Determinatio'i of Diffuser Efficiency 61 

APPENDIX B - Design of Diffuser Model (2) 63 

APPENDIX C - Computer Program Omega 67 

Cl - Assembly of Input Data 67 

C2 - Output 68 

APPENDIX D - Reduction of Data from Sovran and Klomp, Reference 2 . 75 



3 




♦ 



I 



BIBLIOGRAPHY 



80 



FORM DD 1473 



83 



4 



•h ••V.llftll 




LIST OF TABLES 



Table 




Page 


I 


Diffuser descriptions 


35 


II 


Model (1) and Model (2) test data 


51 


III 


Meridional static pressure values for 
Model (1) and Model (2) 


58 


Cl 


Listing of Program Omega 


70 


C2 


Sample output of Program Omega 


74 


D1 


Listing of reference data reduction program 


76 


D2 


Sovran and Klomp, GMR-511, data 


77 



5 



LIST OF ILLUSTRATIONS 



Figure Page 

1 Example of diffusers which previously had 14 

the same description 

2 One-dimensional flow in a diverging channel 17 

3 T-s diagram for an adiabatic, polytropic 20 

compression in a stationary duct 

4 Schematic of the air delivery system 27 

5 Transition ducting, looking downstream 29 

6 Support assembly, including support struts 29 

7 Diffuser support assembly and converging 30 

outer wall (with Model (1) attached) 

8 Diffuser support assembly with Model (1) 30 

cantilevered center-body installed 

9 Center-body of Model (2) 31 

10 Outer wall of Model (2) with support assembly 31 

11 Schematic of diffuser assembly 32 

12 Schematic of diffuser Model (1) 34 

13 Overall view of the experimental apparatus 36 

14 Diffuser inlet plane, velocity survey in 41 

radial direction (measured by hot-wire 

anemometer) 

15 Diffuser inlet plane, velocity survey in 42 

radial direction (measured by five-hole 

pressure probe) 

16 n vs . n for Reference 2 data 46 

17 Q c^ vs. 0, Model (1) test data 49 

18 Comparison of Reference 2 data with 50 

Model (1) and Model (2) test data 

19 Diffuser polytropic efficiency vs. 54 

Reynolds number 



7 




J' 










20 


Dimensionless pressure distribution along 
outer vail of Model (1) 


55 


21 


Dimensionless pressure distribution along 
the inner and outer vails of Model (2) 


56 


22 


Design of Model (2) 


64 


23 


Fitting of circular arcs 


65 


24 


Input field specifications for Program Omega 


69 



8 




♦ ^ 



LIST OF SYMBOLS 



Symbol 



A 



C 



C 

pr 

C . 
pri 






P 



d 

Sc 

K 

L 

L 

L 

w 

T 

M 

m 

mV 

n 

"t 

P 

Ap 



Definition 

Area 

Wetted perimeter 

Coefficient of pressure recovery 

Ideal coefficient of pressure recovery 

Local coefficient of friction 

Average local coefficient of friction 

Specific heat of air at constant 
pressure 

Distance from outer wall 

Unit conversion factor 

Texas Instrument pressure gauge 
conversion factor 

Axial length 

Meridional length along mean 
streamline 

Meridional wall length 

Meridional length along wall 

Mach number 

Mass flow rate 

Referred mass flow rate 

Temperature reading 

Poly tropic exponent 

Total pressure, absolute 

Static pressure, absolute 

Pressure difference between two 
points in the flow 



Units 

. 2 
in 

, 2 
in 



BTU/lb -°R 
m 

in 

ft-lb /lb-sec^ 
m 

in H2O/ count 

in 

in 

in 

in 

slugs/sec 

millivolts 

Ib/in^ 

Ib/in^ 

Ib/in^ 



9 






iil'. 

1 1 1 








q 



R 

g 

R 

e 

s 



T 

t 



c 




V 

w 



X 

7 

e 

C 

X 

0 

T 

$ 

( 

n 



Dynamic pressure (P^ - p) 

Gas constant 
Reynolds number 
Entropy 

Total temperature 
Static temperature 
Temperature 
Temperature 
Average velocity 
Flow rate 

Equivalent flow rate 
Factor (p^ - p^)/p^ 

Ratio of specific heats 
Density coefficient 
Dimensionless static pressure 
Polytropic efficiency 
Dimensionless length 
Viscosity 
Density 

Wall shear stress 
Flow function 
Shape parameter 



Subscripts 

1 Station at diffuser inlet plane 

Ir Referred to station one 

2 Station at diffuser exit plane 

A Atmospheric 



Ib/in^ 

ft - lb/ lb -°R 



o 

o 

o 

o 



R 

R 

C 

F 



ft/sec 

lb /sec 
m 




lb -sec/ft^ 
m 

3 

slugs/ft 

Ib/ft^ 



10 



d Downstream from nozzle 

H Hub or inner wall 

is Isentropic 

p Polytropic 

T Tip or outer wall 

TI Values measured on the Texas Instrument 

Fused Quartz Pressure Gage 

t Total 

u Upstream from nozzle 

w At or along wall 



11 



ACKNOWLEDGEMENT 



The author wishes to acknowledge the guidance given by Professor 
M. H. Vavra of the Department of Aeronautics of the Naval Postgraduate 
School 5 and assistance of Mr. James E. Hammer and the technical staff 
of the Turbopropulsion Laboratory in the preparation of the experi- 
mental portion of this thesis. 

Special thanks go to Mr. Gordon Gulbranson, non-metal model maker 
of the Department of Aeronautics for the beautiful diffuser models 
which he fabricated. 



12 



I. INTRODUCTION 



In the design of annular diffusers, such as those found in jet 
engines, and turbomachines , there has been a problem defining param- 
eters vhich indicate sufficiently the optimum design range for varying 
geometries. Many schemes have been devised, such as the so-called 
'^equivalent conical angle method", to relate a diffuser to a known 
shape. The difficulty in these approaches has been that each person 
using the method developed a slightly different approach which resulted 
in a variety of descriptions for the same diffuser. Gleason [1] gives 
a summs.ry of such methods and points out the difficulty in correlating 
data from different sources. 

Some apparent problems arise when describing a diffuser by inlet 
area, exit area, length, difference in inlet radii, equivalent cone 
angle 3 or combinations of these parameters developed by empirical means. 
Although these terms give some intuitive feel for the way a diffuser 
would perform, two diffusers which differ widely in performance could 
have the same definition. 

The simple example of Figure 1, where the two diffusers have the 
same inlet area, exit area, length and possibly the same divergence 
angle, makes it immediately obvious that the two diffusers would not 
have the same losses. The velocity and pressure distributions through 
the diffuser channel have a direct effect on the boundary layer 
development, which in turn controls the frictional losses. The 
widely differing geometries of the two above-mentioned diffusers 
would cause large differences in boundary layer growth. Therefore, 
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Example of diffusers which previously had the 
same description. 
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the performance of the two would differ and a similar description 
would not be appropriate. 

Sovran and Klomp [2] have empirically developed a fairly successful 
set of parameters for straight wall annular diffusers. The geometric 
area ratio was taken as the primary factor affecting diffusion, and a 
non-dimens ionalized length when combined with the area ratio, gave an 
overall criterion for the pressure-gradient affecting the boundary 
layer development. Their description of annular diffusers is at 
present limited to designs with straight walls. 

Vavra [3 and 4] presented the idea of a shape factor or parameter 
which includes the effects of diffusion and boundary layer influence 
for the evaluation and preliminary design of diffusers of arbitrary 
shapes . 

The objectives of this study were to present the derivation of the 
shape parameter, Q, to substantiate its usefulness from reference data, 
and to evaluate the validity of the results by construction and investi- 
gation of a family of straight wall diffusers such as those in Ref. 2. 

To show that the design shape parameter may be applied to diffusers 
other than those with straight walls, a curved wall diffuser such as 
might be found at the exit of the high pressure compressor of a jet 
engine was designed for an optimum Q and evaluated by tests. 
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II, DEVELOPMENT OF THE SHAPE PARAMETER OMEGA (Q) 



Tbe function of a diffuser is to convert a portion of the kinetic 
energy of the free stream into static pressure by a compression process 
similar to that in a mechanical compressor. The shape parameter^ 

Omega (Q) , to be used for diffuser design , can be evolved by examining 
the differential entropy changes of such a polytropic compression 
process . 

A. ENTROPY CHANGE OF A ONE-DIMENSIONAL, FRICTIONAL FLOW 

An ^expression for the differential entropy change of a polytropic 
compression may be found by investigating a one-dimensional flow in a 
diverging channel, such as shown in Fig. 2, A momentum balance on the 
differential volume of Fig. 2 gives 

m (V / dV) - mV = pA - (p / dp) (A / dA) / (p / dp/2) dA - jCdL 
which reduces to 

m dV = - Adp - TCdL (2) 

or 

_ in tC _ _ . 

dp = - - dV - — dL (3) 

where the values of A, C, p, dL, V, and T are as defined in Fig. 2. 

Introducing the equation of continuity, 

m = pVA (A) 

into equation (3) results in 

dp = - pVdV - (5) 

or 

f = - d(v'/2) - I f dL (6) 
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FIGURE 2 

One -dimensional flow in a diverging channel. 
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with the First Law of Thermodynamics 




or 



^ = dh - Tds 



( 8 ) 



P 



Equating (6) and (7) gives 



dh / d(V^/2) - Tds = - V dL 

p A 



(9) 



The first two terms of (9) combine to the derivative of the total 
enthalpy, H, and if the flow is restricted to an adiabatic case, with 
no energy in any form transferred between the fluid and its surroundings, 
the total enthalpy is constant and its derivative is equal to zero. 

Therefore , 



and equation (8) simplifies to 

Tds = - 7 dL (10) 

p A 

Introducing the concept of the local coefficient of skin friction 
in a one-dimensional flow. 



equation (10) becomes 



Introducing equation (4) and the equation of state for a perfect gas 
in (12) results in an expression for the differential entropy change. 



dH = d (h / V^/2) = 0 



2 C 

Tds = ~ V 7 dL 
2 A 



( 12 ) 




( 13 ) 
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B. EOTROPY CHANGE FROM THER^IODYNAMIC CONSIDERATIONS 

The differential change in entropy can also be found for an 
adiab-atic, polytropic compression process as shown in the T-s diagram 
of Fig. 3. The definition of polytropic efficiency is, 

dT. 

The isentropic relation for a differential compression from p to 
p + dp is 

2 :J: 

T + dT, ^ , -N 7 

, (ziM) 05 ) 

By expanding (dp/p) by a binomial series, and assuming dp to be small 
compared to p, 

dT. A I ^ (16) 

IS 7 p 

Combining (14) and (16) gives 



dT = 1- ZLl ^ 
rip 7 P 



Solving equation (7) for 



,, dT 1 , 



and introducing 



c = R 
P g 7-1 



and the perfect gas law 

p/p = Rg T 

results in an expression for 



(17) 



ds = R K f 

g L 7 -I T p J 

Substituting (17) into (18) gives 




(18) 



( 19 ) 
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T-s diagram for an adiabatic, polytropic compression in 
a stationary duct. 



20 



I 



i 

1 






# 



C. DEFINITION OF OMEGA (fi) 



Equating the two values of differential entropy change from 
equations (13) and (19) results in 



c . R R 

1 2 _E HI. . -s 



2 P ^3 P 



[^ - l] dp (20) 

Ip 

Introducing the concept of referred flow rate, a dimensionless quantity, 

m (21) 

and substituting (21) into (20) and rearranging gives 



■ ,K!i 

™lr p^A^ 






For a polytropic process 

p/p^ = constant 
where n is the polytropic exponent. 
Rewriting (22) in the form 



T ™lr (l^) a 7 ■ l] (£) ^ 



gives two independent differentials which can be integrated. 
Defining the left side of (23) as dX^ and integrating. 



^ A ^ 

I - r!l • 2 A\ C_ 

‘1 J 2 "'ir y A, 



dL 



(24) 



where L is the meridional length along the mean streamline. 

Assuming an average value of coefficient of friction, c^, for the 
whole surface and that m^^^ is a constant for a given flow rate, (24) may 
be rewritten as 
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(25) 



A. 3 



’‘i = r *ir /(r) 



0 

The integral of (25) is a parameter which contains the flow area ratio, 
the diffuser length, and the wetted surface, and will be called the 
**diffuser shape parameter", defined by 



" / (r) r 
0 ^ 



(26) 



D. RELATIONS BETWEEN Q AND DIFFUSER PERFORMANCE 
Defining the right-side of (23) as dX^ , 

pj/Pi 



r , .1/n , , 



If n and n are assumed constant, 
'P 



■2 = - '] < 17 ^ 0 <-> 



For 



P2 = Pi + ^P 
the binominal expansion 



P2^1/n +1 ^Pi + ^ 

vp V p , y 



“ 0 - 7 ^) 



1/n + 1 



= 1 + 






^ 4 - - 



Assuming Ap is much less than and omitting higher order terms fr 
the binomial expans ion, equation (28) becomes 



om 



X = [i- - l] ^ 

2 hp J p^^ 



(29) 
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Equating and 






^ f • 2 „ 

F ”lr “ 



With the equation of continuity, equation (21) becomes 

i 

V 



R T 



m, = d,V,A, . 
Ir 111 p^A 
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Substituting into equation (30) 









l/2p^V^ 



2 = Q 



For a constant total temperature, equation (18) becomes 

dP, 

ds = - R 



8 P. 



Thus equation (19) becomes 



dP. 

P" 



P 



Integrating (33) 






) = 


r-- 1 


1 In ( 






J ' 



( 30 ) 



(31) 



(32) 



(33) 



(34) 



Using the assumption that P = P + AP and p. = p^ + Ap, where AP 

^2 t L ^ i L 

and Ap are small compared to and p, respectively, and with the 
truncated power series expansion of In (1 + x) = x, (34) becomes 



or 





(35) 
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2 

Since for an incompressible flow, M , is small compared to unity, 



Ik - ' 



AP 



( 36 ) 



Substituting (36) into (32) gives 

AP. 



Q = 



(37) 



1/2 

In Ref. 2, the Pressure Recovery Coefficient is defined as the ratio 
of the actual pressure rise and the maximum attainable one for the 
particula. flow rate, or 

c H (38) 

P’' 1/2 

Defining the Ideal Pressure Recovery Coefficient as that which 
would result for an isentropic compression producing the theoretical 



pressure p« of Fig. 3, 
th 






P""^ 1/2 p^V^' 



(39) 



For isentropic flow Ap^ = 0, and assuming density to be constant and 
introducing the equation of continuity. 



pri 

From equation (40), 



p/2 - p/2 

p/2 



= 1 



-d-^) 



(40) 



Then 



p/2 - p/2 [1 . Cp^.l 



Pt -Pj - Pi + P/2 - pj - p/2 11-C il 



or 
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(41) 



- P_ 

' P2 , ^ 

p/2 p/2 

Using the definition of C , (38), and equation (37) there is 



nc 



f 



= c . 

pri 



C 

pr 



(42) 



or 




C . - C 

n 



(43) ■ 



Hence, after having previously determined Q, the average frictional 
loss coefficient, c^, can be calculated from the measured quantities 

Pi’ ^1’ ^1’ P2’ ^2‘ 
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III. TEST INSTALLATION AND INSTRUMENTATION 



The apparatus used in this investigation was located in Building 
230 at the Turbopropulsion Laboratories of the Department of Aeronautics, 
Naval Postgraduate School, Monterey, California. 

A. AIR DELIVERY SYSTEM 

A schematic drawing of the installation is shown in Fig. 4. The 
air source was a 200 H.P.', Carrier centrifugal compressor with a maxi- 
mum flow rate of about four pounds mass of air per second at a pressure 
ratio of about two. Gross throttling was accomplished by a main by-pass 
valve ahead of the air coolers. Fine adjustments of the flow rate were 
made by bleeding air from plenum (1) through a secondary by-pass valve. 

The air from plenum (1) discharged through a flow straightener into 
a four inch pipe containing a flow measuring nozzle. The nozzle used 
was the same as that described by Kelly [5], which had a diameter of 
3.022 inches and a diameter ratio, 3, of 0.734. 

The flow was diffused into plenum (2) by a ten degree conical 
diffuser extending downward from the inlet. The discharge of plenum 
(2) was into a rectangular wooden duct containing flow straighteners 
and a total temperature, iron-constantan thermocouple probe. 

The mean stream velocity in the duct at maximum flow rate was 
approximately 20 to 25 feet per second. Transition pieces shown in 
Fig. 5 were used to change the geometry of the flow channel from 
rectangular to circular and to accelerate the flow to the diffuser 
inlet velocity and pressure. 
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Schematic of the air delivery system. 



B. DESCRIPTION OF DIFFUSER MODELS AND ASSEMBLY 



The diffuser assemblies of the two models constructed and their 
supporting assembly are shown in the photographs of Figures 6 through 
10. The supporting assembly consisted of the duct transition pieces 
with a cylindrical portion, shown in Fig. 6, containing three equally 
spaced support struts for a cantilevered center-body, so that no 
obstructions or supports would be necessary in the diffuser flow 
channel. The flow was accelerated to the diffuser entrance in an 
annular space created by the nose of the center-body assembly and a 
converging outer wall. 

A straight section, where the annular area remained constant, was 
provided just downstream of the inlet to provide room for flow measure- 
ments and surveys, as shown in the schematic of Fig. 11. Smooth 
transitions to the diffuser contour were arranged as shown in Figures 
8 and 9. All computations of the diffuser shape parameter included the 
flow area from the inlet plane to the exit plane of the diffuser, 
inclusive of the cylindrical inlet section. 

The diffuser inlet plane was provided with three supports for 
probes, equally spaced around the outer periphery and bisecting the 
angles formed by the supporting struts of the center-body. These 
supports could accommodate both hot-wire and pressure probes for 
radial surveys. One support was used for a Kiel probe for the measuring 
of total pressure at the inlet plane. Three static pressure taps were 
provided in the same plane, 120 degrees apart on the outer wall. Three 
sets of static taps were placed axially downstream from the inlet static 
taps, at two-inch intervals on the outer wall of the first model, and 
on both the inner and outer walls at one-inch intervals on the second 
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FIGURE 5 

Transition ducting, looking downstream. 




FIGURE 6 



Support assembly, including support struts. 
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FIGURE 7 




Diffuser support assembly and 
converging outer wall (with 
Model (1) attached). 



FIGURE 8 

Diffuser support assembly with 
Model (1) cantilevered center-body 
installed . 
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FIGURE 9 

Center-body of Model (2) 




FIGURE 10 



Outer wall of Model (2) with Support Assembly. 
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Schematic of diffuser assembly. 



model. Static pressure readings were found to be equal at all three 
stations around the circumference, and the three taps at each axial 
location were therefore connected together. 

Model (1) had straight walls as illustrated in Fig. 12, with a wall 
divergence angle of 15 degrees for both the inner and outer contours. 
The construction was such that a series of four diffusers with differ- 
ent shape parameters could be produced by simply cutting off the end of 
the model. Model (2) had a more complicated shape, and was designed 
as described in Appendix B to show the validity of the Cl parameter for 
curved-wall diffusers. A summary of the geometrical characteristics 
of both models is given in Table I. 

C, INSTRUMENTATION 

Pressure readings for the flow nozzle and the diffuser inlet were 
made with a Texas Instrument Fused Quartz Pressure Gage, which is 
shown in the overall view of the experimental apparatus of Fig. 13. 

A constant pressure source calibration of the Pressure Gage, to a 
Meriam Vernier Manometer filled with water, gave a constant value of 
0.800 inches of water per count, with a reading accuracy of + 0.001 
counts. This result was identical to that of Beck. [6]. 

Pressure readings for five-hole probe measurements were made on a 
235 cm, eight tube water manometer bank, which can be read to an 
accuracy of + 0.1 cm. Static pressure measurements along the diffuser 
length were made on a 92 inch, 20 tube water manometer board with a 
reading accuracy of + 0.05 inch. All pressure taps were connected to 
their respective reading devices with plastic tubing. A standard 
mercury wall-barometer v;as used for atmospheric pressure readings. 
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Schematic of diffuser Model (1) 



Diffuser Descriptions 
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Section II for definition of Q 
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Overall view of the experimental apparatus. 
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Temperature readings from the thermocouples upstream of the flow 
nozzle, and in the wooden duct, were obtained with a Leeds and Northrup 
portable potentiometer with compensated cold- junct ion , giving a read- 
ing accuracy of + 0.01 mV. Adiabatic flow was assumed from the area 
of temperature measurement to that of pressure measurement. 

Hot-wire measurements were conducted with a two channel, constant 
temperature, hot-wire anemometer manufactured by Applied Science, 
Carmel, California. The wires were 0.00015 inch diameter tungsten with 
a 0.1 inch gap. Surveys were taken radially from the outer to the 
inner wall. 
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IV. EXPERIMENTAL PROCEHURES AND COMPUTATIONS 



All reduction of data was carried out by the Fortran IV computer 
program *’Omega” described in Appendix C. 



A, PRESSURE 

Pressures obtained from the Texas Instrument Fused Quartz Pressure 
Gage were reduced to pounds per square inch by 

p = 0.03613 K p^j (44) 

where K = 0.8, the conversion factor from counts to inches of water, 
p^^ is the pressure read in counts on the Pressure Gage, and the 
constant value, 0.03613, is the conversion factor for inches of water 
to pounds per square inch from Ref. 7. 

Barometric pressure was converted from inches of mercury to pounds 
per square inch by equation (3-2) of Ref. 6. 

0.0001634 t 



p. 1^1 - ^ ^ 0.0001818 t J 



(45) 



where the coefficients of t^, the room temperature in degrees centi- 
grade, are the thermal expansion coefficients of mercury and brass. 
The factor 0.4912 is the conversion factor for inches of mercury at 
O^C to pounds per square inch from Ref. 7. 



B. TEMPERATURE 

Thermocouple readings, recorded in millivolts, were converted to 
degrees Fahrenheit by the temperature relationship, 

t^ = 32.144 / 35.77 mV - 0.4518 mV (46) 
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obtained from a least squares fit of data from Ref. 8. This relation 
holds for temperatures below 100 degrees Fahrenheit, with the cold- 
junction maintained at 32 degrees, or compensated to 32 degrees by the 
standard cell of the Leeds and Northrup measuring apparatus. 



C. FLOW RATE 



After setting the desired flow rate, the fluid temperatures were 
allowed to stabilize before data were taken. Upstream static pressure, 
p^, differential pressure across the flow nozzle, (p^-p^), and the 
upstream total temperature, necessary to compute flow rate 

as described in Appendix A of Ref. 5. 



w 





(47) 



and 

w = 11.30 e f — 

t y p 
’ u 

vt 2 

where w is the equivalent flow rate, which has dimension in . 

The so-called expansion coefficient, e, is 



(A8) 



e 



2.433 



(1-x) 



1.428 



P 



- (1-x) 



0.286 



1 - .305(l-x) 



1.428 



:] 



(49) 



where 




(50) 



and 

= 1 / 0.0015 (51) 

is a correction factor that takes account of the thermal expansion of 
the nozzle throat diameter. 
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D. CHECK OF NOZZLE FLOW RATE WITH SURVEY DATA 



Preliminary surveys at the inlet plane with a hot-wire probe 
indicated that the velocity profiles at the three stations surveyed 
were similar. Figure 14 is a least-squares fit of the data points from 
all three survey stations. The data were taken at varying flow rates 
at each station, and the velocities were found to vary within + 1 
per cent of the mean obtained with the least-square fit. Hence, the 
flow was very nearly axisymmetr ic , and the non-dimensional profiles 
remained unchanged at the different flow rates investigated. The 
boundary layer was very thin at the inlet. The actual thicknesses were 
0.025 inch on the inner wall and 0.75 inch on the outer wall according 
to the measurements. 

To verify the pressure and temperature measurements, a comparison 
was made of the flow obtained from the survey data at the inlet and 
that measured by the flow nozzle. A survey with a five-hole probe was 
made at the inlet plane to obtain the actual values of the velocities 
which are plotted in Fig. 15. Data for the five-hole probe were cor- 
rected by use of the calibration curves supplied by the manufacturer. 
United Sensor and Control Corporation. The velocity profile obtained 
was essentially the same as that of the hot-wire surveys. The hot-wire 
values were used for the boundary layer regions, where the immersion 
corrections for the pressure probe were unreliable. The yaw angles of 
the flow were less than + 1/2^, and the flow had pitch angles between 
three and six degrees, with the velocities pointing outward from the 
center. The greatest velocities occurred at the outer wall. 
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Diffuser inlet plane, velocity survey in radial direction (measured by hot-wire anemometer). 



I 



I 



I 

I 

f 

* 




I 



# 






1 






o 




o 

o 

to 



ftrj 

lij-j 



vj CVJ — 



Ct-J" 

UJ-J 

z<^ 

?3t> 



V 



0 > 

d 



00 

d 



N 

O 



iD 

d 



iD 

d 



o 



lO 

d 



CSJ 

6 



(S3H0N1) “11VA\ d31fK) 30NV1S10‘ p 



42 



Diffuser inlet plane, velocity survey in radial direction (measured by five-hold pressure probe). 



t? 





The five-hole probe survey showed the total pressures and the 
velocities to be within two per cent of the mean values in the turbu- 
lent central region outside the boundary layers. An average Mach 
number of the flow in this central region was determined from the 
isentropic re lat ion ship 

Zli 1/2 



M 



1 





(52) 



The average static temperature at the inlet was determined by 
assuming adiabatic flow from the upstream total temperature probe to 
the inlet plane, which gave the isentropic relationship 



= 1 + ^ 



M, 



(53) 



Hence , 



T, = 



tl 



1 + 



2zl 

2 



M, 



(54) 



Static pressure was found to be nearly uniform across the inlet 
plane, so that density was almost constant at this station. From the 
perfect gas relationship. 



Pi = 



Pl(144) 



1 R g T, 
g c 1 



(55) 



witli the units of slugs per cubic foot. 

From the survey, the velocities were known as a function of radius, 
V(r). Thus, a flow rate was calculated at the inlet plane by numerical 
integration of 



^1 



Wj = 2 n Dj^ J V(r) r dr 



(56) 



"hi 
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assuming axisymmetric flow and no blockage from support strut wakes. 
This flow rate was found to be three to four per cent higher than the 
flow rate obtained from the flow nozzle measurements. This comparison 
was considered satisfactory for the range of measuring capability. 



E . REYNOLDS NUMBER 

The characteristic length for the Reynolds number (Re) was taken 



as the hydraulic radius of the inlet annulus, AR^, as in Ref. 2, or 



AR 



1 



2 (flow area)^ 
(wetted perimeter) ^ 



(57) 



Therefore , 



Re = 



''i 

u 



and, substituting 



V = Pi Ai . Pi Vi n (R^^- 4^) 



(58) 



(59) 



the Reynolds number (Re) can be expressed in dimensionally correct units 
as 



Re 



V (12) 

g ^ n (R^ + 




(60) 



-7 2 

where viscosity, |j, has a representative value of 3.9 x 10 Ib^-sec/ft 
for the operating range of 85 to 95 degrees Fahrenheit. The flow rate, 
w, is in Ib^/sec, and the radii are expressed in inches. 
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V. ANALYSIS OF TEST DATA FROM REFERENCE 2 



Sovran and Klomp [2] made a study of more than one hundred annular 
diffuser configurations, and their published data was of a form suit- 
able to be analyzed for the present purpose. Of the many other works 
published, none furnished a sufficiently wide cross-section of data 
which could be readily applied. Data from Ref. 2 were reduced by the 
Fortran IV program shown in Table Dl. The results are listed in Table 
D2 and plotted on Fig. 16 for a Reynolds number of approximately 6 x 10 
Each set of data points represents a family of straight wall diffusers 
having the same inlet area and wall divergence angles as shown in 
Fig. 16. The data points fall within a band on the plot of Q c^ versus 
n, which is bounded by lines representing a dimensionless length 
description of the diffuser, L/lSR. An investigation of the values of 
L/AR for existing diffuser designs indicates that the boundaries shown 
on Fig. 16 are in the range of practical designs. The points falling 
below L/Z\R =1.5 represent designs which become so short that diffusion 
is degraded, and those above L/ AR = 13 become excessively long for 
practical use. It can be seen, that as Q falls below approximately 
five, there is an abrupt rise in the value of Q c^, which is repre- 
sentative of the entropy rise and thus of an increase of the losses in 
the diffuser. The trend of the curve for larger values of Q was found 
by considering that, as a diffuser with small divergence angles becomes 
long, its characteristics approach t.iat of a pipe. Reference 9 gives 
the value of the coefficient of friction for an annular pipe as 
c^ = 0.01 at Re = 6 X 10^. A pipe with L/AR = 13 has a value of 
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(PRODUCT OF THE SHAPE PARAMETER AND COEFFICIENT FRICTION) 




FIGURE 16 

n c j vs . n for Reference 2 data. 
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Q = 26 and Q = 0.26, which is a higher value of Q than that 
occurring from diffuser tests at Q = 11.5 and L/ AR = 13, suggesting 
an upward trend in the curve for high values of Q. This trend 
indicates that a minimum of Q c^ occurs between about Cl - 5 and 
Q = 15 , possibly at a value of Q = 10, for an L/ AR = 13. 

Vavra [3] calculated values of Q c^, for two-dimensional subsonic 
diffusers, which lie within the same band as the data plotted on Fig. 
16. 
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VI. DISCUSSION OF TEST DATA 



The test results obtained from diffuser Model (1) and Model (2) are 
given in Table II. The quantities Q for Model (1) are plotted on 
Fig. 17 versus Q. Tests were made at different flow rates to investi- 
gate the effect of Reynolds number on the location of points in the 
n vs. n plot for the different configurations of Model (1). Figure 
17 shows that as Reynolds number increases the quantity Q c^ decreases, 
and vice versa. This effect could be anticipated from the trend of c^ 
against Re for flows in annular pipes, as shown in Ref. 9. 

The curve for L/AR = 6 and the curve for the family of diffusers 
with = 15°, from Fig. 16, and the curves from Fig. 17 were 

transposed to Fig. 18. The shape of the curves from data of Model (1) 
approximates that of the similar family of diffusers from Fig, 16 for 

$ = $. = 15°. However, the values obtained from the tests of Model (1) 

o 1 

are considerably lower than those from Ref. 2. This discrepancy may be 
explained by the fact that the models of Ref. 2 were constructed of wood 
and, therefore, had a higher wall roughness than the cast resin models 
used in this study. The value of c^ for a wood surface would be 
larger than that for a resin surface with lower wall roughness, similar 
to what can be seen from the Moody diagram for pipes of differing 
roughness , 

The values of Q for configuration four of Model (1) became 

slightly negative. Values of fi c^ from Ref. 2 data were also negative 

when C became greater than C .. For an adiabatic process, the value 
pr ^ pri ^ 

n c- cannot be negative, since C is less than C . by definition. An 
f " pr pri 

investigation of the pressure measurements for configuration four of 
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FIGURE 17 

flOf ■ vs. Jl, MODEL (I) test data 



Curve from Sovran S KIomp 
Oota, Rg 16 
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Comparison of Reference 2 data with Model (1) and Model (2) test data. 



TABLE II 



Model (1) and Model (2) test data 



MODEL 


CONFIGURATION 


n 


Q 


% 


R X 10^ 

e 


1 


1 


II. 521 


0.0614 


.8793 


1.165 








0.0667 


.8837 


.879 








0.0566 


.8709 


1.486 








0.0632 


.8628 


1.514 








0.0573 


.8696 


1.518 








0.0569 


.8699 


1.510 








0.0610 


.8794 


1.201 








0.0560 


.8704 


1.526 








0.0859 


.8647 


.768 




2 


9.689 


0.0338 


.8677 


1.517 








0.0356 


.8683 


1.506 








0.0418 


.8734 


1.234 








0.0486 


.8808 


.775 








0.0529 


.8740 


.762 




3 


8.012 


0.0288 


.8699 


.766 








0.0135 


.8633 


1.509 








0.0170 


.8703 


1.208 








0.0200 


.8780 


.980 








0.0146 


.8587 


1.493 








0.0137 


.8607 


1.495 




4 


6.470 


-0.0162 


.8589 


1.502 








-0.0027 


.8756 


.755 








-0.0105 


.8772 


1.034 








-0.0128 


.8677 


1.265 








-0.0151 


.8589 


1.494 
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TABLE II, continued 



MODEL 


CONFIGURATION 


n 


n 


'n 

p 


R X 10^ 
e 


2 


1 


5.701 


0.0949 


.8226 


1.467 








0.1009 


.8279 


1.143 








0.1036 


.8362 


.741 








0.0941 


.8243 


1.466 
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Model (1) showed a fluctuation of the instrument readings of sufficient 

magnitude for C to be larger than C 

pr ^ pri 

It can be anticipated that the curves of Fig. 16 will be shifted 
vertically up or down with Reynolds number similar to the condition in 
Fig. 17. The values of Q for the point of sharply increasing n 
would not be affected by a change in Re, since the entire band of data 
points would be displaced in the vertical direction only. To avoid 
large values of fi a minimum value of = 5 should be used for most 
practical designs. 

To demonstrate the usefulness of the design parameter n> for com- 
plicated shapes, Model (2) was designed as outlined in Appendix B for 
a value of Cl = 5,7, and L/ AR = 6.85. Values of Q = 6 and L/ AR = 6 
were chosen for the original design point, to obtain an n = 0.1. 

By extrapolating from the curve for L/ AR = 6 of Fig. 16 to the 
data of Fig. 17 to show the effect of Re and the difference in the 
model surfaces, the heavy dashed curve of Fig. 18 was obtained. The 
data of Model (2), plotted on Fig. 18, verified the location of the 
L/AR curve for an arbitrary shaped diffuser. 

The efficiencies for both Model (1) and Model (2), plotted on Fig. 
19, remained essentially constant through the Re range of this study. 

Static pressure (p^) on the inner and outer walls of Model (1) and 
Model (2) was non-dimens ionalized by 




and plotted against non-dimensionalized , meridional wall length, 

X = oo Figures 20 and 21, where 1 is the meridional distance 

along tile wall from the inlet to the pressure tap, and is the 
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DIFFUSER POLYTROPIC EFFICIENCY VS. Re 
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X (Dimensionless Wail Length) 

FIGURE 21 
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meridional wall length. Values for X, and are contained in 
Table III, for the outer wall of Model (1) and for both the inner and 
outer walls of Model (2). The static pressure became higher than the 
ambient pressure at about X = 0.5 on the inner wall of Model (2), 
which gives a Q less than that at the exit. The value of Q reaches a 
minimum between X = 0.7 and X=0.9 and increases again toward Q = 0.293 
at the exit plane. This condition may be due to a separation of the 
flow, and indicates the need for proper wall contouring in addition to 
satisfying the criterion of Q. 
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TABLE III 



Meridional, static pressure values for Model (1) and Model (2) 



CONTOUR 


STATION 


X* 


p (in. water) 
gage 


c** 


OUTER 


1 


0.0 


-13.54 


1.0 


CONFIG. 1 
Model (1) 


2 


0.133 


- 8.79 


0.729 




3 


0.267 


- 5.46 


0.540 




4 


0.400 


- 3.50 


0.429 




5 


0.533 


- 2.24 


0.357 




6 


0.667 


- 1.45 


0.312 




7 


0.800 


- 0.82 


0.276 




8 


0.933 


- 0.34 


0.249 




9 


1.000 


0.0 


0.230 


OUTER 


1 


0.0 


- 9.66 


1.0 


CONFIG. 4 
Model (1) 


2 


0.4 


- 4.90 


0.720 




3 


0.8 


- 1.45 


0.520 




4 


1.0 


0.0 


0.435 



* X - The dimensionless length 1/L 

w 



**C “ Dimensionless pressure 
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TABLE III, continued 



CONTOUR 


STATION 


X 


p (in. water) 
w 


c 


INNER 


2 


0.107 


-7.97 


0.787 


Model (2) 


3 


0.222 


-4.27 


0.558 




4 


0.336 


-2.14 


0.426 




5 


0.460 


-0.61 


0.331 




6 


0.578 


+0.21 


0.280 




7 


0.697 


+0.73 


0.248 




8 


0.817 


+1.01 


0.231 




9 


0.933 


+0.85 


0.241 




10 


1.000 


0.00 


0.293 


OUTER 


1 


0.0 


-11.41 


1.000 


Model (2) 


2 


0.121 


-10.91 


0.969 




3 


0.241 


-7.90 


0.783 




4 


0.366 


-5.18 


0.614 




5 


0.484 


-3.50 


0.510 




6 


0.605 


-2.40 


0.442 




7 


0.718 


-1.58 


0.391 




8 


0.823 


-1.15 


0.364 




9 


0.938 


-0.37 


0.316 




10 


1.000 


0.00 


0.293 
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VII. CONCLUSIONS AND RECOMMENDATIONS 



This paper presents a design parameter Q> for diffuser designs, 
which takes into account losses due to the geometry of the diffuser. 

It is shown from reference data and model tests that there is an 
abrupt increase in losses for diffusers designed for a value of Cl 
less than about five, and that this effect is independent of Reynolds 
number . 

The parameter Q is defined for incompressible, adiabatic flows. It 
is recommended that further experiments be made to extend the data for 
incompressible flows to higher values of Cl, and that an investigation 
of the compressible case be made. 

It is suggested that Model (2) be improved by recontouring the 
inner wall to avoid separation. It is recommended also that the 
present inner contour be cut off at the station where the value of Q 
is a minimum. Addition of static pressure taps on the inner wail at 
the inlet plane and a method for surveying the exit plane would pro- 
vide data to better describe the flow. 

It is recommended that the method of measuring pressure be investi- 
gated and improved to give more precise results for small, fluctuating 
pressure measurements, prior to conducting further tests on diffusers 
with Q c^ approaching zero. 
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APPENDIX A 



DETERMINATION OF DIFFUSER EFFICIENCY 



The diffuser efficiency was determined by the method of Vavra [10] 
for an adiabatic compression process with friction. The polytropic 
efficiency is defined for the differential compression from p to 
p / Ap of Figure 3 as 



dT. 

_ IS 



(Al) 



The dimensionless flow function and the relation for the poly- 
tropic exponent, n, from Ref. 10 are 

1 



R 



1 



$ - w 1—^ T . 

C (Sc 1 hh V 



2/n 



n-f-l 



22_ 






1 



and 



n = 



1-7 (1-T1p) 



(A3) 



where w, T^ , Tj^, and p^ are either measured or computed for the 



particular diffuser design. 

For an incremental change in pressure from p^ to p^^ the pressure 
P 2 may be expressed as p^ / Ap, where Ap is small with respect to p^. 
The pressure ratios in (A2) can be expressed as 



Pi V p,y 



(A4) 



By the binominal expansion theorem, 

2/n 2/n 

0 



n p- 



(A5) 
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and similarly 



n-fl 
rPox ^ 



(:2y ^ ... 

Vpj^y n 



(A6) 



and by assuming Ap to be small, the (Ap/p^) and higher order terms are 



much smaller than Ap/p^^ and are omitted. 



From (A2) 



2/n 



n+1 



c 27 



vp^y Vp^y 



(A7) 



Substituting expressions (A5) and (A6), and solving for n, gives 

T 



^ [2 - i1 

Pi L Tj J 

n — - rn 



pi |2 . S + I + -^ 

27 c rp p. 



(A8) 



Equating (A3) and (A8) gives 



with 






^Gr- 



(A9) 



^ = 1 + ^ M? 



2 1 



the polytropic efficiency is 

2 ^ [1 + (7-1) M^] 






- 7M^ [2 -^ + 1 
c ' 1 L p. 



(A 10) 



(All) 
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APPENDIX B 



DESIGN OF DIFFUSER MODEL (2) 

From the results of Ref. 2 and the test data of Model (1), a design 
point of Q = 6 and L/AR = 6 was chosen from Fig, 18 for Model (2). The 
inlet area of Model (2) was the same as for Model (1). 

From equation (42) of the main text, 

C . - C = n 



pri 



pr 



or, from equation (40), and a value of Q = 0.1, 



[l .(^] . O.l + c 



pr 



(Bl) 



For the range of C from 0.65 to 0.72 for a L//iiR = 6, equation (Bl) 

pr 

gives 






(B2) 



2 1 - (0.1 + C ) 

pr" 

2 

giving a range of from 39 to 51 in . 

An approximate radius was chosen to give a reasonable axial opening 

2 

at a radial exit plane for = 45 in . 

Since the outer surface, E, of Fig. 22 was most critical from the 
standpoint of flow separation, it was designed first. A smooth curve 
of the desired shape was sketched in for surface E. To approximate 
the curve sketched, a series of circular arcs were constructed which 
provide for simple tooling layouts. 

Figure 23 shows the method by which two circular arcs are fitted to 
satisfy the condition that both have centers on the same line at their 
tangent point (4), and that they fit into a given length, L. 
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FIGURE 22 

Design of Model (?) 
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FIGURE 23 

Fitting of circular arcs 
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The largest radius, R^, is chosen to give the desired effect of 
the original sketch. A curve is then drawn from point (1) to point 
(2), where (2) is an approximate point of the desired curvature change. 
Point (2) defines a and b for any given distance d. The center of R 2 
becomes point (3). 

From geometry, R 2 and R *2 can be defined as 



and 




(B3) 




Point (2) may be iterated until R 2 = ^^ 2 * centers of R^ and R 2 

will then be on the same radius at their tangent point and the curve 
will fall within the axial length L and height d. 

The same procedure as described above was extended to surfaces of 
more than two circular arc segments, such as surface E, of Fig. 19, by 
locating the arc centers E-1, E-2, and E-3 , Surface F was constructed 
in a similar manner. 

With areas and known, a linear distribution of flow area from 
station (1) to station (2) was chosen. Circles of the diameter which 
gave the desired area at a given meridional station were then con- 
structed along an approximate mean streamline as shown by the example 
circles on Fig. 22. Surface B was then sketched in, tangent to the 
area circles, and fitted with circular arcs as described above. The 
flow channel was then checked for uniformly increasing area as a function 
of length, and adjusted to prevent any convergence in the flow area. 
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APPENDIX C 



COMPUTER PROGRAM OMEGA 

Program "Omega^^ is a modification of program "Diffuser” used by 
Beck [6]. The program is written in Fortran IV for the IBM 360-67 
digital computer. The program was established to compute the values 
of flow rate 5 Re, Q, and r\^ by the methods described in the main 

text. Table Cl is a listing of the program. 

Cl. Assembly of Input Data 

The field specifications for the input data are shown in Fig. 24. 
Card one of the data deck is a 72 space "Hollerith" statement for data 
identification. The content of this card is printed immediately below 
the output title. Card two is a two digit I-format number indicating 
the number of data points in the deck. Card three is the input of 
inlet and exit diameters and the value of Cl. Cards four and five con- 
tain the data recorded for one data point and must be repeated for 
each point. No other cards need be repeated. 

All floating point numbers may be placed anywhere in their field. 
I-format numbers must be right justified. 

The definitions of the program input variables are given below. 
Fortran 



Symbol 


Definition 










Format 


NUM 


Number of sets- of 


data 


points 




I 


2 


DOl 


Diameter of 


outer 


wall 


station 


one 


F 


10.0 


Dll 


Diameter of 


inner 


wall 


station 


one 


F 


10.0 


D02 


Diameter of 


outer 


wall 


station 


two 


F 


10.0 
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Fortran 

Symbol 


Definition 


Format 


D02 


Diameter of outer wall, station two 


F 


10.0 


DI2 


Diameter of inner wall, station two 


F 


10.0 


OMEGA 


Value of shape parameter 


F 


10.0 


IRUN 


Number of run 


I 


3 


MO 


Month 


I 


3 


IDAY 


Day 


I 


3 


PATM 


Atmospheric pressure 


F 


10.0 


PINOZ 


Static pressure upstream of nozzle 


F 


10.0 


DPNOZ 


Pressure differential across nozzle 


F 


10.0 


TNOZ 


Total temperature upstream of nozzle 


F 


10.0 


TEMP 


Temperature read at barometer 


F 


10.0 


PI 


Static pressure at station one 


F 


10.0 


Q1 


Dynamic head at station one 


F 


10.0 


TTl 


Total temperature measured upstream 
of station one 


F 


10.0 


TRM 


Ambient room temperature at exit 


F 


10.0 


C2. Output 









The program output is in three sections. The first section is a 
printout of data as read from the input cards as a check. The second 
section is the calculation of flow rate which is divided into input 
values converted by the program to dimensionally correct units, and 
output of flow rates computed by the program. The third section is 
also divided into input values converted by the program and output 
values of interest in evaluation of diffuser performance. 

Table C2 is a sample output of the program. 
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CARD 



1 72 




72 Space Hollerith Statement 




Data Identification 


Blank 



CARD 2 



1 2 




I 3 




NUM 


Blank 



CARD 3 



1 10 


11 20 


21 30 


31 40 


41 50 




FIO.O 


FIO.O 


FIO.O 


FIO.O 


FIO.O 




DOl 


Dll 


D02 


DI2 


OMEGA 


Blank 



CARD 4 



1 3 


4 6 


7 9 


10 19 


20 29 
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APPENDIX D 



REDUCTION OF DATA FROM SOVRAN AND KLOIP, REFERENCE 2 
Data of Ref. 2 were reduced by a Fortran IV program to a form 
useable for this study. Table Dl is a listing of the program. The 
diffuser description of Ref. 2 was the same as in Fig. 10. 

The definition of the input variables are 



Fortran Symbol 


Description 


AD?AH 


The divergence angle, 0^, of the inner wall 


ALFAT 


The divergence angle, o^, of the outer wall 


DIMLL 


The dimensionless length, L/AR. 


ARATIO 


Area ratio, A^/A^^ 


RRATIO 


Radius ratio, 


RTl 


Radius of outer wall, station one 


CPR 


Coefficient of pressure recovery 



The axial length was computed from the geometry and a numerical 
integration was made to obtain Q. The ideal pressure recovery 
coefficient was computed from the area ratio and used to compute Q 
The results are contained in Table D2 , where 



Fortran Symbol 


Description 


OIEGA 


The shape parameter, Q 


CFBAR 


Average local coefficient of friction, c^ 


DIMLL 


The dimensionless length, L/AR 


OCF 


The product. Cl c^ 
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Listing of reference data reduction program 
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TABLE D 2 



SOVRAN & KLQMPt GNR-511, DATA 



OMEGA 


CFBAR 


3o 54 5A 


Co0027 


5o2663 


0o0054 


7o2076 


0o0117 


9ol634 


OoOT 25 


9o9462 


0o0136 


10o6331 


O 0 OI 37 


llo52r7 


Co 013 2 


3o6454 


C 0 OOI 3 


7o2076 


O 0 OII 9 


7o2C^6 


Co Cl 24 


7o2C76 


O 0 OI 43 


9o 1634 


O 0 OI 35 


9ol634 


O 0 CI 43 


llo5207 


O 0 OI 47 


2o 96 39 


O 0 O 22 I 


3o4l99 


O 0 O 345 


3o7999 


0o0291 


4o4838 


Co036<= 


4o5598 


0o0291 


4o 6358 


0o0406 


5o0918 


O 0 O 39 I 


5o3198 


0o0372 


5o5478 


0oC388 


2o9639 


CoC262 


3o7999 


0o0?06 


3o7999 


0o0391 


3o7999 


O 0 O 375 


3o 7999 


0o0388 


3o7999 


0o0430 


4o5598 


0o0296 


4o5598 


O 0 O 344 


4o5598 


0o0368 


4o5598 


0o0421 


2o71^2 


O 0 O 5 II 


3ol363 


O 0 O 575 


3o4570 


0o0631 


2o42^0 


0oC957 


2o7^30 


O 0 IO 57 


2o98ll 


0ol098 


3o4800 


Co 1070 


3o7C44 


0ol046 


3o7914 


0o0983 


2o7884 


O 0 O 5 I 8 


2o9829 


O 0 O 547 


3o0478 


0o0570 


3o3071 


0o0551 


3o5017 


0o0556 


3o6314 


0o0514 


3o8908 


0o0471 


4o 1501 


0o0427 



DIMLL 


OCF 


2,2800 


0,0097 


3o 6300 


0,0283 


5o6600 


0,0341 


80 3700 


0,1146 


9,7200 


0,1351 


lloC709 


0, 1461 


13,1000 


0,1519 


2o28C0 


0,0047 


5o6600 


0,0861 


5, 6600 


0,0891 


5o5600 


Col031 


8,3700 


0,1236 


80 3 700 


0,130 6 


13,1000 


C, 168<= 


2o260C 


0,0655 


2o93C0 


0,1181 


3o6100 


0,1105 


4o960C 


0,1 656 


5o6400 


0,1328 


6o3200 


0,188^ 


80 3 400 


0,1993 


10,3700 


0,1978 


13,CB00 


0,2151 


2o260C 


0,0775 


3,6100 


0,1165 


3,6100 


0,1435 


3o6\C0 


0,1425 


3,6100 


0,1475 


3,6100 


Oo 1635 


5,6400 


0,1348 


5,6400 


0,1568 


5,6400 


0,1678 


5,6400 


0,1918 


2,2500 


0,1387 


2,9300 


0,1804 


3,6000 


0,218? 


2,2400 


0,2 320 


2,9100 


0, 2900 


3,5900 


0,3274 


6,2900 


0,3722 


9,6800 


0,3875 


13,0600 


0,3728 


2,3600 


0,1444 


2,7100 


0,1633 


3,0500 


0,1736 


3,4000 


0,1856 


3,7500 


0 , 1 949 


4,4400 


0,1866 


5,1400 


0,1833 


5,8300 


0,1772 
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TABLE D2 , continued 



SOVRAN & KLOMP, GMR-511, DATA 



OMEGA 


CFBAR 


DIMLL 


OCF 


2o7884 


0o0582 


2o?600 


Oo 162^ 


2o9829 


0o0624 


2o7100 


0ol863 


2o2947 


Ool ^08 


2o3400 


0o323i 


2o4?58 


0ol429 


2, 6800 


0o34R0 


2o5578 


0ol468 


3o0300 


0o3756 


2o6602 


0ol48^ 


3o38C0 


Oo 39-^9 


2o7467 


Col515 


3o7300 


0o4l63 


2 o 88 20 


0ol48/- 


4o420C 


0o4?76 


2o9835 


0ol^43 


5ol2C0 


0o43C5 


3o05S0 


0ol397 


5o800C 


0o4273 


2ol696 


Col691 


2o3300 


Co3669 


2o2974 


0ol723 


2 o 6800 


Co3959 


2o4030 


0ol738 


3o0300 


0o4177 


2o48e6 


0ol761 


3o3700 


0o4382 


2o5627 


0ol768 


3o7200 


0o4530 


lo974B 


Oo218e 


2o32C0 


Co4321 


2o0752 


0o2247 


2 o 6600 


Co4662 


2ol5C3 


0o2278 


3o01C0 


0o4919 


2 o 2281 


0o2309 


3o3600 


Oo 5146 


2o2851 


Oo 2?22 


3o 7100 


Co 5305 


2o4784 


0o2310 


5o7900 


0o5725 


2o85O0 


0o0378 


2o4200 


O 0 IO 8 I 


3o8475 


0o0298 


3o9200 


Oo 1 1 47 


4o7895 


0o0249 


60 IIOO 


0oll90 


5o5813 


O 0 O 2 IO 


9o 1700 


O 0 II 72 


lo8581 


-0oCC30 


lo4300 


-OoOOS'^ 


2o45 53 


O 0 OI 53 


2o3000 


0o0376 


3o0217 


0o0276 


3o6000 


0o0833 


3 o 47 34 


0a0328 


5o3400 


C 0 II 39 


lo8581 


-OoOC36 


lo4300 


-Co0C67 


lo8581 


O 0 OI 26 


lo4300 


O 0 O 233 


2o4553 


O 0 OI 24 


2o3000 


O 0 O 3 O 6 


3o0217 


0o0226 


3o6000 


0o0683 


3o4734 


0oC230 


5o34CC 


O 0 O 799 


3o4734 


0o0369 


5o3400 


O 0 I 28 O 


2o3698 


-C 0 CI 55 


lo^700 


-0o0368 


3o3697 


-O 0 CO 57 


2o3400 


-0o0191 


4o4930 


-Co 00 12 


3o6^0C 


-O 0 OO 53 


5o5348 


0o0C37 


5o3800 


O 0 O 2 O 6 


2o3698 


-C 0 OIO 9 


lo4700 


-0o0258 


5o5348 


O 0 OIO 2 


5o3800 


0o0566 


lo6569 


-O 0 OOO 4 


lo500C 


-O 0 OOO 6 


lo 8824 


OoOOOl 


lo9500 


O 0 OOO 2 


2oC503 


O 0 OOO 3 


2o4000 


O 0 OOO 6 


2o1776 


OoCOOO 


2o8500 


OoOOOO 


2o2757 


-O 0 OOO 2 


3o3000 


-O 0 OOO 4 


2o3525 


OoOCOl 


3o 7500 


O 0 OOC 2 


2o4134 


O 0 I 735 


4o2000 


Co4186 


2o4631 


Oo 1^23 


^0 660 0 


0o4244 


2o5024 


0ol725 


SollOO 


Co4316 
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TABLE D2 , continued 



SOVRAN & KLOMP, GMR-5U, DATA 



OMEGA 


CFBAR 


DIMLL 


OCF 


2o53A6 


O 0 I 732 


5o5600 


Co 4390 


lo6569 


0ol409 


lo5000 


0,233^ 


2o5?A6 


0ol848 


5o56C0 


0o4680 


2o8i.54 


O 0 CO 2 I 


2 oC 100 


OoOC58 


3o5666 


O 0 OO 45 


2 o 8600 


0o0162 


4o2207 


O 0 OIO 8 


3o8100 


C 0 O 457 


4o7081 


O 0 OI 22 


4o7100 


0o057^ 


5o4215 


O 0 OI 35 


60 5 200 


C.,0734 


5o6828 


O 0 OI 35 


7o^200 


0o0768 


5o 8977 


OoO] 20 


8o310C 


0o0705 


2o81 54 


O 0 OO 33 


2 o 0100 


O 0 OIC 8 


2o8053 


0o0067 


2 . 0 OOOO 


C 0 OI 88 


lo4279 


Oo2348 


lo560C 


Oo3353 


lo5757 


0o3059 


2 o 0200 


Co4820 


I 068 O 3 


0o3143 


2o4900 


0o5281 


lo7532 


0o3205 


2o9500 


0o5619 


I 08 O 75 


Oo3232 


3o4200 


0o5841 


lo 8472 


0o3242 


3o6800 


0o5989 


lo4279 


Co2313 


lo5600 


Oo 3303 


1o8472 


Oo3291 


3,8800 


C 06 O 79 


2o2C^7 


O 0 II 29 


2oC800 


0o2488 


2o4275 


0o0844 


2o5400 


0o2048 


2o857o 


0o0903 


3,9300 


0o2582 


2o2C47 


O 0 I 2 OI 


2 o 0 800 


0o2648 


2o8579 


OolU3 


3o 9300 


0o3182 
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